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. The influence of the electric field created by a gate potential of the silicon 
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q . quantum computer on the hyperfine interaction constant (HIC) is obtained. 
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Q , The errors due to technological inaccuracy of location of donor atoms under 

a gate are evaluated. The energy spectra of electron- nuclear spin system 
of two interacting donor atoms with various values of HIC are calculated. 
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The presence of two pairs of anticrossing levels in the ground electronic state 
is shown. Parameters of the structure at which errors rate can be greatly 
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minimized are found. 
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C : I. INTRODUCTION 

cr. 

The first approach to the NMR quantum computer realization simultaneously proposed in 
1997 by two groups of investigators and then experimentally confirmed @,|], consists 
in using organic liquids, where different molecules with a few connected with each other 
nonequivalent nuclear spins-qubits function as independent quantum computers. Another 
approach, which is qualitatively different from the bulk-ensemble one and is still unrealized, 
was proposed and analyzed in details by Kane || and elaborated in [§,13,14]- It is based on 



creating artificial multi-spin system and using individual addressing to different spins-qubits. 
For this purpose it is suggested to use silicon-based structure of MOS-type, where donor 
atoms of stable phosphorus isotope 31 P are implanted into the thin layer of spinfree silicon 
isotope 28 Si at a definite depth. Donor atoms replace silicon ones in knots of crystal lattice. 
Such a donor has shallow impurity states, which have a large magnitude of the effective Bohr 
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radius and nuclear spin I = 1/2. Every donor atom with nuclear spin in semiconductor 
structure is supposed to be placed regularly with adequate accuracy under its "own" control 
metallic gate (gate A), which is separated from the silicon surface by a thin dielectric layer 
(for example, silicon oxide with the thickness of about some nanometers). A-gates form 
linear regular structure of arbitrarily length with period 1 (Fig. 1). Changing of the electric 
potential of J-gates, which are placed between A-gates, allows us by redistributing electron 
density between neighboring donors to control the degree of overlapping of electron wave 
functions, which are localized on neighboring donors a and b, and to control the constant 
of exchanging interaction J, and the constant of scalar interaction of their nuclear spins 
I a b as well. It is assumed that with the help of the electric field, induced by A-gates, one 
can change the distribution of electronic density near the nuclei in the ground state by 
choosing individually, correspondingly, each donor atom nuclear spin resonance frequency, 
which is determined by the hyperfine interaction with its electron spin. This allows to fulfill 
quantum operations by the way of selective reaction of resonance radio-frequency pulses on 
nuclear spins of the given donors. The second section of the work studies the influence of 
the electric field, created by potential V on the A-gate, on hyperfine interaction constant 
(HIC) of the donor atom in two cases: for gates having the form of a round disc and that 
of a strip. In the third section HIC fault approximation is made, which is determined by 
the technological placement inaccuracy of donor atom under the gate. In the fourth section 
the results of full calculation of the electron-nuclear spin system energy spectrum of two 
interacting donor atoms with different magnitude of HIC are demonstrated. Separately the 
phenomenon of hyperfine energy levels crossing and anti-crossing is investigated, owing to 
which the information about nuclear spin state is transferred to the electron spin of donor 
atom when moving through the point of crossing. This property was suggested to be used 
in the process of electrical measurement of a definite nuclear spin state • It is also very 



important when using magnetically-mechanic resonance method [12 
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II. DEPENDENCE OF THE HYPERFINE INTERACTION CONSTANT ON THE 

ELECTRICAL FIELD 



Spin Hamiltonian of the hyperfine interaction is: 

H IS = A (It) , A (V) = y |* (0, V)\ 2 • (£) (1) 

where /Ujv = 5.05 ■ 10~ 27 J/T — nuclear magneton, /i = 47r ■ 1CT 1 T 2 sm 3 /J, g N = 2,26 
- Lande's factor for 31 P, fig = 9.27 ■ 10~ 24 J/T — Bohr magneton, and ^ Q (0, 1 7 ) — wave 
function of the electron in the ground state placed on the nuclei. Experiments give the value 
of the probability density of the electron to be on the nuclei (0)) 2 = 0.43 • 10 24 sm" 3 at 

v=o i. 

Taking into account that intrinsic semiconductor at low temperatures behaves like di- 
electric and supposing that D >> I a ~ c» d (Fig. 1), we express the electric field on the 
donor atom depending on potential V on gate A, (round disc with radius a — -^), which lies 
on the surface of the dielectric. To find the electric potential induced by the round disc at 
the point with the coordinates r* = (p', z') , p' 2 = x' 2 + y' 2 , we use known expression 0: 



/ / A 2V 

(p[p ,z ) = — arctg 

IT 



2a 2 

(2) 
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\ p' 2 + z' 2 -a 2 + yV 2 + z' 2 - a 2 ) 2 + Aa 2 z 

On the line, which passes through the gate center p' = \Jx 12 + y' 2 = , near the donor 
atom (z' = c), we receive: 

<p (p', z>) = —arctg- - E c (z' - c ) + ^{z>- cf - ^-p' 2 + ..., (3) 

71 C I I 

where 

2V a W ac _ AV ^ oc 4 

C vr ' a 2 + c 2 ' ^ tt ' (a 2 + c 2 ) 2 ' c vr V ( Q 2 + C 2)V2 

the electric field, its gradient in the z' axis direction and the electric field gradient along 
the radial direction near the donor atom. When crystalline lattice and external electric field 
are absent, the wave function of this state is trans- formed into the hydrogen-like function, 
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which corresponds to IS state. Shroedinger equation for this function in the electric field 
must be complemented by the perturbation operator AH — e ■ tp (p', z') ( e = 1.6 ■ 10~ 19 C 
- elementary charge) ||. Using ([|), (|]) in x, y, z basis (Fig. 1) we receive the following 
expression (p 2 = x 2 + y 2 ): 

?Vp n pF/ pR" . . . . 

AH = —arctg- + eE c ■ z - —± ■ z 2 + — ^ • p 2 = AH + AH X + AH 2 + AH 3 (5) 

7T c 2 2 

The wave function Fj (r) correction, specified by the perturbance (|5|), is: 

Fj = if s) + Aif s) + Aif 2 s) , (6) 

where 



A _ AH 2s ,ls (2s) (7 s 



- the first order correction of the perturbation theory, AH 2s i s = J d 3 r ■ 
(Fj 2s ^ (r)) ^AH^j Fj ls ^ (r) - - matrix element of the perturbation operator and SE = 
Ei s — E 2s — energy residual of unperturbed donor states, and 

lAff I 2 

M,2 S) =-VS 7i f\* > m = 2p > 3p ' • • • (8) 

- the second order correction (here E m — unperturbed energy levels of donor, AH m ^ n 
- matrix elements of the perturbation operator). During further calculations we shall use 



hydrogen-like functions as unperturbed wave functions. They are independent of j [TO 



(r) = 1 —n,e-^; (r) = — =J m ( 2 - - ) e"* (9) 

v^F- (a*f 2 Ay/tot- (a*f 2 \ 1 
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where r 2 = x 2 + y 2 + z 2 (Fig. 1), and a* = 4 ^:"a = 20 A — effective Bohr radius 
@, e = 11.9, e = 8.85 • 1(T 14 F/sm - - vacuum permittivity, h = 6.62 • 1(T 34 J*sec, 
m* = 0.31m e , m e = 9.1 ■ 1CT 31 kg. The second order terms in give an expression a the 
respective correction of the hyperfine interaction constant of the donor atom: 



AA (V) _ (F) 2 - (F^) ^ AF W AF 2 (ls) ( AF[ 



(18)' 
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The second term in (fLOf), according to |8|Jl^|, may be evaluated as: 



^V = --" 4^o (a*) 3 (11) 
where (AE) 1 — some mean value of (E m — Ei s ) 1 . At AE ~ 0.04 ev, c = 2a = 10 

2AF l ' ls ' 

nm it was obtained that (V in volts) F(1 2 3) = — 0.19V . 

When calculating Aif 2 s,is , we shall take into account that the item in the perturbance 
AHq does not contribute because of orthogonality of functions F^ (r) and F^ (r) , and 
AHi does not contribute because of oddness of the function by z under integral. Resulting, 



AH 2s , ls = ^eE' c (a*) 2 - (a*) 2 (12) 



Taking into account that |^ = ^ , (|T0|) appears as: 



— 4— = - 9TCo (a } ~ae ¥5E { 1 -¥¥ C ) + {¥^E—{ 1 ~¥¥ C )) (13) 

The energy residual 5E can be evaluated, if we use the formula for hydrogen-like atoms 
8E = — | • , e2 ; . The evaluation for silicon gives SE = —0.023 ev. Using the above- 
mentioned evaluations with the help of ( |i~3"D for c = 2a = 10 nm we shall result: 



AACV) „ 
M = 0.55V -0.09V 2 



-- --■ < 14 > 

Let us consider the case, when the A-gates have a form of infinitely long strips with the 
width 2a = Ia ■ They are placed at the distance D from the conducting substrate, from 
which the gate potential is counted out. Conducting similiar operations we receive for the 
donor atom HIC the following expression (c = 2a = 10 nm, — = 100 ): 

- -0.063V 2 (15) 



Voltage in (15) is expressed in volts, as written above. 
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III. INFLUENCE OF THE TECHNOLOGICAL SPREAD OF DONOR ATOMS 
ON THE HIC. VOLTAGE ERROR ON A GATE 



Technological inaccuracies of the donor atom place under the gate relatively x = 0, z = c 
are designated Sx and Sz . We investigate only the strip gate variant: Factorizing potential 
expression for strip gate by degrees of 5x and 5z , we write it as: 

f (*, z) « v (0, c) + *" (Sx g + fc) (, - c) + &g + fc > (z - c)> + 3lf ^± fc) ^ 

(16) 



where 



9y>(fa,c+fe) _ _ p J, _ (fa) 2 (2c 2 -a 2 ) 



gVffe^c+fe) _ p, L fe(2 C 2-a 2 ) (fa) 2 (4c 4 +a 2 c 2 ~a 4 ) 1 , > 

dz 2 -^cj 1 C (a 2 + C 2) 2 C 2(a 2 + c 2)2 h U ' j 

M^££±M - 1 F" fl fe(2c 2 -a 2 ) (fa) 2 (2e 2 +a 2 ) 1 
d(x 2 ) - 2^c Y c ( a 2+ c 2) 2(a 2 + c 2)2 J 



We consider 5z ~ (5a;) 2 , that is why corrections of (5z) 2 degree are omitted being of the 
next order. Using expressions fll3|) , (|H)[)-(|l7l), we receive the final expression for the relative 
error of the HIC: 

5A . f „ 2c 1 , p n2 f 9 2c 2 -a 2 2c 4 -a 4 ] 

— = 0.063V 2 — \ + (5x ^ 0.063V 2 ? - 0.085V 5 f 18 

A I a 2 + c 2 J V ; \ (a 2 + c 2 ) 2 2c 2 (a 2 + c 2 ) 2 J V ; 

The second item in (|IBD nullifies in a variety of values of a, c and V. That is why by 
defining necessary relative error it is possible to determine structure parameters a and c and 
voltage V, when the needed accuracy is realized. One percent relative error of the hyperfine 
interaction constant corresponds to the donor atom deviation along the z axis, which is 
connected with the technological inaccuracy of placing in 2-3 nm. 

Voltage error on the gate, which is admissible during correctly conducted quantum cal- 
culation process, can be determined from the final expressions for HIC ([14]) and (|l^). The 



strip width of radio frequency pulses must be at least one order smaller than the resonant 
frequency of nuclear spins, which is about hundreds of kHz. Considering the HIC value, 
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which corresponds to the strip width of resonant pulse S (AA) ~ 10 4 Hz, for the gate volt- 
age error we receive 10~ 4 -j- 10 _3 V. It is worth of noting that in the case of round discs the 
distinguished voltage exists, that leads to the significant increasing of the admissible voltage 
error. It is caused by the presence of the linear by V term in (|Hp . 

IV. ELECTRON-NUCLEUS SPIN SYSTEM ENERGY SPECTRUM OF TWO 

INTERACTING DONOR ATOMS 

Consider two donor atoms separated from each other by distance 1. This distance must 
satisfy the following condition: constant J of the electrons' effective exchange interaction, 
determined by the partial overlap of their wave functions, must correspond to the possibility 
of making two-qubit operations. Spin Hamiltonian of this system is: 



H = 2fi B B ^S a +S b j+J ^S a S b j +AH = H + AH, 
AH = -g N fi N B (f a + + A a (f a s)j + A b (f b S b "j , 

where S a , I a , S b , I b - spin operators of the electron and the nuclei for the first and the 
second atom, B — the magnetic field (it is directed in parallel to z axis (see Fig. 1)), A a 
and A b — hyperfine interaction constants depending, generally speaking, on gate potentials, 
Hb — Bohr magneton, fi N — nuclear magneton, g N = 2.26 — Lande factor for 31 P . 

We consider \M a ; M b ; m a ; m b ) to be the basis (eigen states in the strongest fields), where 
M a , M b , m a , m b — quantum numbers of the electrons' and nuclei's spin projections. 
Designating them by arrows up (f ) and down (j) , we have (numbers are strings or rows in 
the matrix): 



ll> = 


ITTTT) 


|5> = 


ITITT) 


|9> = 


= IITTT) 


|13> = 


= lUTt) 


|2> = 


ITTTI) 


|6> = 


inn) 


|10) 


= um> 


|14) 


= iun> 


|3> = 


ITTIT) 


|7> = 


ITUT) 


|H) 


Hint) 


|15) 


= unr) 


|4) = 


ITTIi) 


|8> = 


inn> 


|12) 


= UTU) 


|16) 


= 1111!) 



(20) 
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When nuclear spins are absent (AH = ), Hamiltonian ^ has four times four degener- 
ated eigen states (Fig. 2). Two underlying electron levels have the point of crossing C when 
= 1. 

Symmetric matrixl6*16, corresponding to the matrix of eigen states of Hamiltonian 
Ah = ^j- , in the considered basis has only diagonal elements —ji + aa ^ ab , aa ~ ab , — aa ^ ab , 

a _ "a+Qft _i_ Qq-Qh ag+Q-b a a +a b a a -a b a a —a b a a +a b a a +a b < a a —a b 

h 1 4 5 H 1 ' 4 5 4 ' 4 ' 4 5 p 4' 4' 4 ' ~r 4 J 

-/i- 2 ^, - a f st , 2a^ ) ^i+aaja, where // = a a = ^, a fe = ^ and non-diagonal 

elements, from which we point out only up-diagonal elements: 

( A ' l ) 5 .= ( A/! ) H = ( A '0 13 a„=( A/i ) 



15,12 2 



fA^ =fA^ = (Ah) = (Ah) 

V /9,3 V /10,4 V /13,7 V / 



(21) 

_ Qo 

14,8 2 



Since the projection of full electron and nuclear spin in the direction of the magnetic field 
M + m = M a + Mb + m a + m& is conserved, the matrixl6*16 breaks down into 5 reduced 
matrixes of smaller dimension, corresponding to the values M + m = 0, ±1,±2 , which 
diagonalize independently. The next states correspond to these 5 matrixes: 

m + M = : |4),|6),|7),|10),|11),|13) 
m + M = 1 : |2),|3),|5),|9) 

m + M = -1 : |8),|12),|14),|15) (22) 
m + M = 2 : |1) 
m + M = -2 : |16) . 

With the help of the found matrix, introducing Hamiltonian (|i~9D, the picture of the 
hyperfine splitting of two donor atoms' electrons' levels was achieved. It is shown for a a = 
0.3; «b = 0.4 in figure 3. 

The interaction of electronic and nuclear spins not only takes down the four times de- 
generacy, but brings near point C the phenomenon of anticrossing of the levels, which have 
the same value of m + M. 

For two crossing electronic levels it takes place for the states with m + M = — 1 |15 > 
and 1 12 > (two solid bold lines) and with m + M = |10 > and |13 > (two dashed lines) 
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(see Fig. 3). 

As the result of this anticrossing the state |15 > with M a + M b = — 1 and m a + m b = 0, 
when J increases (with the help of J-gate), makes a transition to the \12 > state with 
ma + mb = — 1 and M a + M b = 0, i.e. the simultaneous inversion of one nuclear spin and one 
electron spin takes place (the information about nuclear spin condition is transfered to the 
electron spin). Herewith the electron triplet state (S = 1) makes a transition to the singlet 
one (S = 0). 

The energy residual of the underlying anticrossing state |15 > and crossing state |14 > 
in the second order of perturbation theory for the strong magnetic fields (/3 » 1 ) and for 
A a = A b = A , according to 0, is 

where v j — transition frequency, which corresponds to projection changing of two nuclear 
spins. It increases, when (3 decreases. 

In the levels |15 > and |14 > are suggested to be used as "working" ones to conduct 
the calculation process and to measure the final state. 

However, we also can propose an alternative variant. It is possible to use the second pair 
of anticrossing levels with m + M = 0. Underlying state of the anticrossing level 1 13 > in 
the strong fields corresponds to M a + Mf, = — 1, m a + m& = 1, and above- lying state |14 > 
corresponds to M a + M b = — 1, m a + m b = 0, i.e. they differ by the inverted value of one of 
the nuclear spins. 

When J increases, the state 1 13 >, moving through the anticrossing point, makes a 
transition to the state 1 10 > with M a + M b = and m a + m b = 0, i.e. the simultaneous 
inversion of the electron spin such as the nuclear one happens (the information about nuclear 
spin is transferred to the electron spin, too). 

According to , the measurement process of nuclear spin is s upposed to be implemented 
in two stages. By changing the value of the electron exchange interaction J from J < 2^bB 
up to J > 2fisB the adiabatic transition of the nuclear spin condition is made from triplet 
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electron ground state to the singlet one. Then the electron state measurement is made. 

If at first the electrons of two donors are at the triplet state, simultaneously being 
placed on its own atom, then, when J increases, it is more energy- advantageous for the both 
electrons to exist in the singlet state being placed together on the one of the donors, if the 
bonding energy of two electrons on the single donor exceeds the bonding energy for each 
electron on its own donor (D~ -state). As a result the charge carrying takes place from one 
donor to another. As it is pointed out in H, this fact can be registered by highly sensitive 
single-electron capacity methods. 

At the present time different methods of the spin states measurement are discussed. In 



11] it is suggested to use a scheme, agreeably called "turnstile", which omits current of 



electrons only with the definite spin state. This electric method allows to measure electron 



spin states. Another approach was discussed in |12| . It is based on using magnetic resonant 
force microscopy. 



V. CONCLUSION 

In this work an influence of the electric field, inducing by the A-gate (for the gate as a 
round disc and as a strap) on the hyperfine interaction constant was analyzed in details. 
The weaker dependence of the hyperfine interaction constant from the depth of underlying 
donors under the gate, and also decreasing of its magnitude when the distance from the gate 
to the substrate increases are noted. Conducted calculations with two gate types allows 
us to make a conclusion about significant dependence of the hyperfine interaction constant 
from the form of a gate. 

Parameters a and c (5-10 nm) of the structure at which errors rate can be minimized 
are determined by (17). Obtained expressions allow us to determine admissible values of 
technological inaccuracies 5x and 5z (2-3 nm) when the constant A errors and working 
potential values V (0.1-1 volts) of a gate are defined. In work the full energy spectrum of the 
electron-nuclear spin system of two interacting donor atoms for A a ^ A b was calculated. The 
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existence of two pairs of anticrossing levels at the ground electron state, corresponding to the 
full electron and nuclear spin projection m a +m b + M a + M b = —1 and m a +m b + M a + M b = 
was demonstrated. Underlying components of these pairs correspond to the transition from 
the triplet electron state (S=l) to the singlet one (S=0), when the field decreases. If the 
full electron spin projection changes from M a + M b = — 1 to M a + M b = 0, the nuclear spin 
projection also correspondingly changes. Thus the information about nuclear spin state 
transfers to the electron system. Both pairs of anticrossing levels may be used to conduct 
calculation and measurement process. Herewith we must use the frequency, determined by 
splitting between underlying anticrossing levels (|13 > or |15 >) and a bove-lying level |14 >, 
to change nuclear spin state. 

In conclusion authors would like to thank L.A. Openov for the valuable comments. 
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